OiJOt,FV KNOX LIBRARY 

" vATE SCHOOL 
WON a3&43-5101 



ll<: 



m 







I 




i 




Approved for public release; distribution is unlimited. 



Three-Dimensional Effects of Crack Closure in Laminated 
Composite Plates Subjected to Bending Loads 

by 

Mehmet Baskaya 

Lieutenant Junior Grade, Turkish Navy 
B.S., Turkish Naval Academy, 1988 

Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE IN MECHANICAL ENGINEERING 

from the 

NAVAL POSTGRADUATE SCHOOL 
June 1994 



11 



REPORT DOCUMENTATION PAGE 



Form Approved OMB No. 07()4-()l«K 



Public rep^»ninj; burden lor thu» cullecUon of inlormaUon is esumaied lo average 1 hour per response, including the lime for reviewing instnicijon. searching existing data 
sources, gaihcnng and maintaining the data needed, and completing and reviewing the collection of infonnatjon Send comments regarding this burden estimate or any 



d.cr i 1 thCi vxlivCtion ct ud^tuiaUon, tricluJmg sdggcaiivtiS fcf leducibg this burdeij, Washington Scfvicts. Oirscl^alt (or InlormaUoii OpcTaticms 

and Reports, 1215 Jefferson Oavis Highway, Suite 1204, Arlmglon. VA 22202-4302, and to the Office of ManagemeDl and Budget Paperwork Reduction Project (0704- 
01 XX) Washington DC 20503 



1 A(}ENCY USE ONLY 


2. REPORT DATE 


3. REPORT TYPE AND DATES COVERED 




June 1994 


Master’s Thesis 



4 TITLE AND SUBTITLE THREE-DIMENSIONAL EFFECTS OF 
CRACK CLOSURE IN LAMINATED COMPOSITE PLATES 
SUBJECTED TO BENDING LOADS (UNCLASSIFIED) 


5. FUNDING NUMBERS 


6 AUTHOR(S) Mehmet Baskaya 


7 PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 
Naval Postgraduate School 
Monterey CA 93943-5000 


8. PERFORMING 
ORGANIZATION 
REPORT NUMBER 


9 SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


1 1. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the 
official policy or position of the Department of Defense or the U.S. Government. 


12a DISTRIBimON/AVAILABlUTY STATEMENT 

Approved for public release; distribution is unlimited. 


12b. DISTRIBUTION CODE 
*A 



1 ,v abstract 



Fracture is one of the dominant failure modes in structures subjected to external loads. Stress and 
deformation fields around the crack tip are important to understand the crack propogation and arrest. For a plate 
w ith a through-the-thickness crack and subjected to a bending load, there is crack closure on the compression side 
of the crack face. The present study investigates effects of crack closure on the stress and deformation fields on 
the tension side of the crack face. A three-dimensional finite element analysis is performed for laminated 
composite plates using both the line and surface crack closure models. For a composite whose longitudinal elastic 
; modulus is much greater than the transverse modulus, line and surface closure models result in higher stresses 
j near the crack tip in comparison to the no-closure solution. Hence, no-closure solutions are nonconservative for 
the composite. Transverse shear is the major cause for the nonconservative solution. 



14 SUBJECT TERMS FINITE ELEMENT METHOD, BOUNDARY 
CONDITIONS, CRACK CLOSURE MODELING 



17 SECURITY 

CLASSIFICATION 
OF REPORT 



18. SECURITY 

CLASSIFICATION OF 
THIS PAGE 



19. SECURITY 

CLASSIFICATION 
OF ABSTRACT 



15. NUMBER OF 
PAGES 4 6 

16. PRICE CODE 

20. LIMITATION 
OF ABSTRACT 
UL 



Unclassified 



Unclassified 



Unclassified 



NSN 7540-01-280-5500 



Standard Form 298 (Rev. 2-89) 

Prescribed by ANSI Sid. 239-18 



I 



ABSTRACT 



Fracture is one of the dominant failure modes in 
structures subjected to external loads. Stress and deformation 
fields around the crack tip are important to understand the 
crack propogation and arrest. For a plate with a through- the- 
thickness crack and subjected to a bending load, there is 
crack closure on the compression side of the crack face. The 
present study investigates effects of crack closure on the 
stress and deformation fields on the tension side of the crack 
face. A three-dimensional finite element analysis is performed 
for laminated composite plates using both the line and surface 
crack closure models. For a composite whose longitudinal 
elastic modulus is much greater than the transverse modulus, 
line and surface closure models result in higher stresses near 
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I . INTRODUCTION 

Plates and shells are very commonly used as structural 
members in construction of ships and other structures. Due 
to their wide usage as structural members, much research has 
been performed to understand their failure mechanisms and 
modes under different loading conditions. Fracture is one of 
the most frequent failure modes in plates and shells. 

The first analysis on plates containing centrally 
located through- the- thickness cracks and subjected to 
bending loads was performed by Williams [Ref. 1] in 1951. In 
his analysis, Williams used the forth-order classical plate 
bending theory with Kirchhoff boundary conditions. He showed 
that there are some discrepencies between the stress 
distributions for extensional and bending problems. In 
William's bending study, surface shear stresses on crack 
faces didn't vanish. This physical violation is attributed 
to the fact that the Kirchoff bending theory can not satisfy 
the physical boundary conditions on the crack surface 
completely. 

This problem was solved by Knowles and Wang [Ref. 2] . 
They used Reissner/Mindlin ' s sixth order plate theory which 
allows three boundary conditions at edges. They solved the 
problem by employing a singular integral equation. However, 
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their solutions were limited to plates with very small 
thicknesses. Analyses including effects of the plate 
thickness on the bending stress distribution were performed 
independently by Hartranft and Sih [Ref. 3] and Wang [Ref. 
4] . 

All the theories mentioned above assume that the crack 
surface is stress -free and straight. In reality, for a plate 
under a bending load, there exists crack closure on the 
compression side and crack opening on the tension side of 
the crack surface. A photoelastic experimental study 
performed by Smith and Smith [Ref. 5] showed partial crack 
closure. A significant change of local stresses was observed 
as a result of the crack closure. 

A finite element analysis was first undertaken by Jones 
[Ref. 6] using Kirchhoff's bending theory to study effects 
of crack closure in plate bending. He compared the closure 
and no- closure solutions, and reported increases in crack 
opening displacements and intensity factors due to existing 
crack closure. In Heming's study [Ref. 7], the 
Mindlin/Reissner sixth order plate theory was employed, 
which is more realistic and allows satisfaction of three 
boundary conditions at the crack face. These results were 
contradictory to the results obtained by Jones [Ref. 6]. 
Heming stated that there were significant decreases in the 
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bending stresses and displacements on the tension side due 
crack closure. 

In those studies, the line closure model, which assumed 
the outermost points from the neutral axis on the 
compressive side of the crack surface to be in contact, was 
used. However, Smith and Smith's experiment [Ref. 5] showed 
that there was closure over an area of the crack surface, 
the size of which varied depending on the plate and crack 
dimensions. A three dimensional finite element analysis was 
performed by Alwar and Nambissan [Ref. 8] to investigate the 
closure region and effects on the stress distribution of 
this closure. A very crude finite element mesh was used in 
their analyses. 

Finally Kwon [Ref. 9] studied the crack closure effect 
for unidirectional composite plates using the line closure 
model and the Mindlin/Reissner plate theory. He found that, 
the crack closure on the compression side of a bent 
unidirectional composite plate increased the crack opening 
displacement and bending stress on the tensile side. This 
finding was in contradiction to that for an isotropic plate. 
Hence, it was concluded that results obtained without 
considering crack closure were nonconservative for 
unidirectional composite plates. The nonconservative 
solutions for unidirectional composite plates were caused by 
the effect of transverse shear. 
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In the present study, effects of surface closure on 
stresses and displacements around the crack tip in laminated 
composite plates are investigated using three-dimensional 
finite element analyses. 
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II. FINITE ELEMENT FORIOTLATION 



A. OVERVIEW 

A finite element program is developed to model a plate 
with centrally located through- the- thickness crack and 
subjected to a constant bending load. As shown in Figure 1, 
the cartesian coordinate system is employed in the model. 
Displacements in the x, y and z directions are represented 
by u, V, and w respectively. 




Figure 1 Geometric representation of the plate model 
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The place is assumed Co have simply supporced boundary 
condicions along ics edges. Due Co symmecries along Che x 
and z axes, only a guarcer of Che place is analyzed Co 
reduce Che number of nodes needed in modeling Che place. 
Therefore, Cwo symmeCric boundary condicions and Cwo simply 
supporced boundary condicions are required for Che quarcer 
place. Three-dimensional parallel -piped elemencs are used in 
conscruccing Che sCiffness macrix. Each node has Chree 
degrees of freedom and, Cherefore, chere are 24 dofs per 
elemenc as shown in Figure 2 . 




Figure 2 Nodal coordinaces of a Cypical finice elemenc 
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Displacements of the plate are functions of x, y and z 



coordinates. That is 









( 1 ) 



To interpolate the displacements within each finite element, 
shape functions are utilized. Displacement can be expressed 
as , 



where subscript i represents the node number, and the 
numeric subscripts 1, 2, and 3 denote the displacements in 
the X, y, and z directions, respectively, and is a 

function of x, y and z. 

From the elasticity theory, strains are expressed as 
given below. 



8 



8 8 



// = S N, * Un 



v= Z N,*U,-2 Z 

i=\ i=\ 



( 2 ) 



1=1 






( 3 ) 





, dv , dw 
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When written in the matrix form, 



{e}6„ =[^]6.24{«}24.1 (4) 

where {e} is the strain vector of size 6 by 1, {u}® is a 
nodal displacement vector of size 24 by 1, and [B] is a 
matrix containing derivatives of the shape functions. 

The three-dimensional generalized Hooke's law is 
defined as, 

{a}=[D]{e} (5) 

where [D] is a 6 by 6 material property matrix, and {a} and 
{e} are stress and strain vectors. The stress and strain 

vectors are written as, 

{o} = [Oxx, ^yy> Ozz, Xxyi "^xz, "ryz] 

( 6 ) 

{ e } = [e joc » ^yy / e zz / Yxy , Yxz , Yyz] ^ 

The elastic strain energy can be written as, 

£ = i I {£) ’‘[D]{e)dJ'= i j { (/) U}dV (7) 

V V 

In equation (7) , only the terms in parenthesis have 



8 



dependency on Che coordinate system. The volume integral of 
these terms yields the stiffness matrix, 

=\[BV[D][BW ( 8 ) 

r 

The equation of static equilibrium for an element can 
be written as, 

= ( 9 ) 

where {F}® is the vector of forces acting on the elemental 
nodes . 

B. CONSTITUTIVE EQUATION 

The material property matrix [D] is a symmetric matrix 
which includes elasticity moduli and Poisson's ratios of the 
material. In the present work, isotropic, orthotropic and 
laminated composite plates are examined. 

1. ISOTROPIC MATERIAL 

There are three independent material constants in 
the [D] matrix for an isotropic material. Elasticity moduli 
and Poisson's ratios are the same in every direction. For an 
isotropic material matrix [D] is given as. 
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( 10 ) 



£*(l-u) 




E*v 


( l+u)*( l-2*u) 
E*Vi 


(l+u)«(l-2*u) 

£»(l-u) 


( l+u)*( l-2*u) 
£*u 


(l+u)*(l-2*u) 

£-*u 


(l+\j)*(l“2*u) 

E*u 


(l+u)»(l-2»u) 

E*v 


(l+u)»(l-2«u) 


(l+u)»(l-2»u) 


(l+u)«(l-2«u) 


0 


0 


0 


0 


0 


0 


0 


0 


0 



0 0 0 
0 0 0 
0 0 0 

GOO 
0 G 0 
0 0 G 



2. ORTHOTROPIC MATERIAL 

The number of independent material constants 
in the [D] matrix is six for an orthoropic material. The 
material properties have three principal axes . The material 

properties of the plate are considered as, Ej/Ej=5, 20, 40, 

and 100 , Gj,/E,=0.6, Gjj/E,= 0.5, Uj^=0.25, and u, 3 = 1 / 9 . Inverse 

of the [D] matrix is called the compliance matrix and it is 
given in [Ref. 10] as shown below. 



\!E\ -U2i/£2 -U3i/£^3 

-\>nlE\ MEi -U32/£'3 
-Un/fi -023/^2 1/^3 

0 0 0 
0 0 0 
0 0 0 



0 0 0 

0 0 0 

0 0 0 

I/G23 0 0 

0 I/G31 0 

0 0 I/G12 



( 11 ) 



where E denotes inplane elasticity modulus, G denotes shear 
modulus, and niomeric subscripts 1 , 2, and 3 represent 
material principal axes. 
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3. LAMINATED COMPOSITE MATERIAL 

A laminated composite plate with [0/90] ^ is 
investigated. There are two [D] matrices in the model to 
represent the orientations of 0 and 90 degrees. The matrix 
[D] for the 0 degree layer is the same as that for the 90 
degree layer except that material properties in the x and z 
directions (see Figure 1) are interchanged for the 90 degree 
layer. As a result, the compliance matrix for the 90 degree 
layer is given below, 



C. STRAIN-DISPLACEMENT RELATION 

Linear shape functions are utilized in the present 
finite element formulation. That is. 



MEy -U23/£2 -y>n!E\ 0 0 0 

—\)32^E'i X/Ej —\j\2/E\ 0 0 0 

— U3i/£3 —\J 2 \IE 2 ME\ 0 0 0 



( 12 ) 



0 0 0 I/G 21 0 0 

0 0 0 0 l/Gn 0 

0 0 0 0 0 I/G 32 




(I) • (x) • e?) 
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( 13 ) 




/^6 = (l)-(^)-(f) 
Nl = (^) • (f) ♦ (I) 
A'8 = (i) ♦©•(§) 



where a, b, and c are dimensions of a finite element as seen 
in Figure 2 . 

Matrix [B] , which relates strains to nodal 

displacements of an element is formed by taking partial 
derivatives of the shape functions. Hence, six strain 

components are expressed as shown below, 

8 

<«= Z |[A', •£/,!] 

/ = 1 

8 

S i\N,*Va] 

/ = 1 

8 

<■== Z i[A'.»W/3l (H) 

/ = 1 
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8 8 

<-V= I |[M*Wil+ S 

/ = 1 / = 1 

8 8 

/ = 1 / = 1 

8 8 

I |[M*f/,2]+ Z i^[N,*Ua] 

/ = 1 / = 1 

D. STIFFNESS MATRIX 

The element stiffness matrix [K]^ is found by volume 
integration of the term ([B]'^[D] [B] ) as in the equation (7), 

where x varies from 0 to a, ( x € [o a] ) 

y varies from 0 to b, { y E [o b] ) 

z varies from 0 to c, ( z G [o c] ) 

The global stiffness matrix is constructed by adding 
element stiffness matrices. In the present finite element 

mesh, the number of nodes in the x, y, and z directions are 
21, 9, and 21 respectively. It gives the halfbandwidth of 

636 and total degrees of fredom of 11907. A static pressure 
load is applied at the nodes ao the xz plane lying in the 
middle of the plate. 
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III. BOUNDARY CONDITIONS AND MODELING OF CRACK CLOSURE 



Selecting proper constraints imposed along the edges of 
a plate is very important. As mentioned in [Ref. 7], 
application of Kirchoff’s boundary conditions or 
Reisner/Mindlin boundary conditions results in contradictory 
solutions. It is more realistic to employ the 
Reissner/Mindlin theory, which describes three boundary 
conditions along the edges of a plate. Therefore, in the 
present models the Reissner/Mindlin BCs are utilized. 
Because of the symmetry involved, a quarter of the plate is 
modeled. The model has two symmetric and two simply 
supported boundary conditions. Symmetric BCs are satisfied 
by constraining displacements in the normal direction to the 
boundary surface. Simply supported boundary conditions, on 
the other hand are imposed by constraints in vertical 
displacements and tangential displacements at the boundary 
surface . 

Crack closure is a constraint on the crack face. There 
are three crack surface models which are used in analyses of 
plates with through-the-thickness cracks. The first model 
assumes that the plate has no constraint on the crack 
surface, as seen in Figure 3a. This is called the 
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Figure 3a No-Closure Model 




Figure 3b Line Closure Model 




Figure 3c Surface Closure Model 
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no- closure model. The crack surface is considered as 



stress- free and straigth . 

The assumption of no- closure allows crack opening to 
occur on the tension side of the crack face. On the other 
side of the crack face, deformations of the same magnitudes 
occur in the opposite direction. Hence, overlapping or 
penetration takes place by employing the no- closure model. 
However, penetration is not possible in the actual case. In 
reality, there would be contact over a region on the crack 
surface. And this contact is called crack closure. 

We can model the closure region in two different ways. 
In the first model the outermost nodes on the compression 
side with respect to the neutral axis are constrained to 
avoid penetration. This is called the line closure model and 
is shown in Figure 3b. Alternatively, a region which 
contains some of the nodes between the outermost nodes and 
neutral axis on the compression side of the crack face can 
be constrained, as shown in Figure 3c. This second model is 
called the surface closure model. 

Plate bending theories permit only the line closure 
model, because the theories assume that plane sections 
remain plain before and after bending. The experimental 
study by Smith and Smith [Ref. 5] shows that crack closure 
takes place over a portion of the crack face. Therefore, the 
surface closure model is the closest describtion of the 
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crack closure. An iteration process is necessary to model 
surface closure using a contact algorithm. Proper number of 
nodes are constrained to ensure that there is no penetration 
on the crack surface. 

When using the finite element method, it is possible to 
model the surface crack closure using 2-D plate bending 
elements or 3-D solid elements. If 2-D plate bending 
elements are used in modeling, a stack of plate elements are 
used through the plate thickness. As a result, the crack 
face can have a piecewise linear deformation and the closure 
portion becomes a rectangular shape. With 3-D solid 
elements, displacements at the crack face can have more 
general deformation. Hence, a general shape of surface crack 
closure can be modeled. Because of this advantage, we use 
3-D elements in the present study. 
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IV. RESULTS AND DISCUSSIONS 



In the present study, cracked plates made of three 
kinds of materials are examined. These are isotropic, 
orthotropic and three layered composite plates with an 
orientation of [0/90/0] . Four different elasticity moduli 
ratios, El/E2=5, 20, 40, and 100, are considered for 

composite and orthotropic models. Consequently, we analyzed 
nine models with different material properties. Three types 
of crack closure modelings, which are discussed in the 
previous chapter, are employed to model the crack surface. 
These are surface closure, line closure, and no- closure 
models. All three models are compared with each other for 
each cracked plate. 

All the stress and displacement data are computed on 
the plane containing the crack. A three dimensional 
representation of the crack surface for the no- closure model 
is shown in Figure 4. In this figure flat section is the 
uncracked portion of the crack plane. Since for the 
no- closure model there is no constraint on the crack 
surface, displacements with the same magnitudes and in the 
opposite directions are observed on the tension and 
compression sides. In the surface closure model, due to the 
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crack closure, only displacements in the crack opening 
direction are allowed, as shown in Figure 5. 

Effects of closure on the crack opening displacements 
can be explained by considering a compressive force due to 
closure as shown in Figure 6a, that acts on the compression 
side of the crack surface. A statically equivalent force and 
moment, as in Figure 6b, can be thought acting in the middle 
of the plane. While the force increases the crack opening 
displacements on the entire surface of the crack, the moment 
causes decreases on the tension side. For smaller elasticity 
ratios the effect of the moment is greater. Hence, we 
observe decreases in crack opening displacements. For higher 
elasticity ratios, especially in the case of a composite 
plate, the effect of the force is higher. 

When the elasticity ratio is higher, shear deformation 
due to transverse shear becomes more important . As seen in 

Figure 6c the rotation caused by the moment is decreased due 
to the transverse shear effect. Here 0 is the rotation angle 
for the moment without considering transverse shear effect, 
and y is the rotation angle in the opposite direction due to 

the transverse shear. 

In the displacement plots, the horizontal axis 
represents the normalized distance from the crack tip and 
the vertical axis represents crack opening displacements on 
the tension side of the plate. These crack opening 
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TENSION SIDE 





COMPRESSION SIDE 

gure 5a Compressive Force due Co Closure 



TENSION SIDE 




gure 6b SCaCically Equivalenc Force and Moment 




gure 5c Transverse Shear Effect 
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displacements (CODs) are normalized by the maximum 
displacement of the corresponding no closure model . Since 
there is contact on the compression side for closure models, 
only displacements on the tension side are plotted. 

It is observed that the crack opening displacements for 
line and surface closure models are very close to each other 
for all the elasticity ratios. When low elasticity ratios 
are considered, crack closure decreases the crack opening 
displacements on the tension side all the way. This is due 
to the fact that the effect of the moment is higher for low 
elasticity ratios. Since there are more constraints for the 
surface closure model its displacements are smaller. Figure 
7 shows the crack opening displacements on the tension side 
of a laminated, cross-ply composite plate with the 
elasticity ratio El/E2=5. 

From Figure 8 it is seen that, when the elasticity 
ratio, E1/E2, is 20, crack closure reduces the displacements 
near the crack tip. On the other hand, it increases the 
displacements away from the crack tip. When the plots for 
El/E2=5 and El/E2=20 are compared, a slight increase near 
the crack tip is also observed. 

Further inreases in the crack opening displacements are 
obtained for elasticity ratio El/E2=40, as seen in Figure 9. 
Near the crack tip, the displacements for closure cases are 
still lower than that for the no closure case. 
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CRACK OPENING DISPUCEMENTS ON THE TENSION SIDE 




NORM. DIST. FROM THE CRACK TIP 



Figure 7 Crack Opening Displacements (CCDs) on 
side of a laminated, cross-ply composite 
elasticity ratio E./E; = 5. 



the tension 
plate with 
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eiascicicy ratio E./E,=20. 



CRACK OPENING DISPLACEMENTS ON THE TENSION SIDE 




NORM. DIST. FROM THE CRACK TIP 



Figure 9 Crack Opening Displacements (CODs) on the tension 
side of a laminated, cross-ply composite plate with 
elasticity ratio E./E,=40. 



26 



CRACK OPENING DISPUCEMENTS ON THE TENSION SIDE 




Figure 10 
side of 
elascicity 



Crack Opening Displacements (CODs) on 
a laminated, cross -ply composite 
ratio E,/E, = 100. 



tne tension 
elate with 



When El/E2=100, closure increases the crack opening 
displacements all the way, as seen in Figure 10. This is 
because the transverse shear effect is high enough to 
overcome the effect of the moment also near the crack tip. 

Since the closure cases reduce the stresses on the 
compression side more, stress distributions on the tension 
side are different from those on the compression side. 
Therefore, both tension and compression sides are examined. 
In the stress plots, the horizontal axis represents the 
normalized distance from the crack tip upto the corner of 
the plate along the uncracked portion of the crack plane. 
The vertical axis represents stresses that are normalized by 
the applied constant pressure . 

From Figures 11 and 12 it can be seen that for El/E2=5, 
there exist large discrepencies in the stress levels between 
closure and no- closure cases. These differences occur both 
on the tension and the compression sides of the plate. On 
the other hand, the differences in stresses for surface 
closure and line closure cases are very small as seen in the 
crack opening displacements. This is especially true on the 
tension side. 

When higher elasticity ratios are considered, increases 
in the stresses near the crack tip are observed for closure 
cases compared to no- closure case, as shown in Figures 13 to 
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18. When the E1/E2=100, higher stresses are obtained on the 
tension side for closure cases. 



29 




E./E, = 5 . 



MAXIMUM STRESSES ON THE TENSION SIDE 




NORM. DIST. FROM THE CRACK TIP 



Figure 13 
iaminaced, 
E./E,=20 . 



Stress discributior. on the tension side of a 
cross-ply composite plate with elasticity ratio 




Stress distribution on the compression side of a 
, cross-ply composite plate with elasticity ratio 



Figure 14 
laminated 




E / £. = 4 0 . 



1 




E./E. = 40 . 
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E./E, = 100 . 



in M 




3_/E2 = 100 . 
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V. CONCLUSIONS AND RECOMMENDATIONS 



When a laminated composite plate with a through the 
thickness crack has a large ratio between the longitudinal 
and transverse elastic moduli, crack closure increases the 
crack opening displacement and tensile bending stress around 
the crack tip compared to those obtained without considering 
crack closure. There is also a slight difference between 
surface and line closure solutions since the surface closure 
model has more constraints. Both crack closure solutions are 
nonconservative. As a result, crack closure must be 
considered for most composite plates under consideration. 

An experimental study shows that the surface closure 
model is the closest description of crack closure. Since 
similar displacements and stresses are, however, observed 
around the crack tip, line closure can be used to model 
crack closure for computational efficiency. 

In future studies experiments may be performed to 
validate the computer results. Effects of crack closure on 
composite plates containing delamination cracks need to be 
analyzed . 
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